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ABSTRACT 

Holmberg IX X-l is a well-known ultraluminous X-ray source with an X-ray luminosity of 
- lO^ergs -1 . The source has been monitored by the X-ray Telescope of Swift regularly. Since 
2009 April, the source has been in an extended low luminosity state. We utilize the co-added spectra 
taken at different luminosity states to study the spectral behavior of the source. Simple power-law 
and multi-color disk blackbody models can be ruled out. The best overall fits, however, are provided 
by a dual thermal model with a cool blackbody and a warm disk blackbody. This suggests that 
Holmberg IX X-l may be a lOM© black hole accreting at 7 times above the Eddington limit or a 
lOOM© maximally rotating black hole accreting at the Eddington limit, and we are observing both 
the inner regions of the accretion disk and outflows from the compact object. 

Subject headings: accretion, accretion disks - binaries: close - stars: individual: Holmberg IX X-l - 
X-rays: binaries 



1. INTRODUCTION 

Ultraluminous X-ray sources (ULXs) are luminous 
(Lx > 10 39 ergs _1 ) nonnuclear X-ray point-like sources 
in galaxies with apparent X-ray luminosities above the 
Eddington limit for a typical stellar-mass (~ 10M©) 
black hole. The majority of ULXs are believed to be 
accreting objects in binary systems due to their strong 
X-ray flux variability observed on timescales of hours to 
years. Assuming an isotropic X-ray emission, ULX is 
the best candidate of intermediate-mass black hole with 
a mass of - 10 2 - 1O 4 M (Makishima et al. 2000; Miller 
& Colbert 2004). While ULX may represent a missing 
link between stellar-mass black hole and supermassive 
black hole in galactic center, its formation and evolution 
is not well understood. 

The X-ray spectral properties may provide some hints 
about the connection between ULXs and Galactic black 
hole X-ray binaries. In particular, many of the ULXs can 
be modeled with a multi-color disk (MCD; Mitsuda et al. 
1984) plus power-law model which is a popular spectral 
model for Galactic black hole binaries (see Remillard & 
McClintock 2006). In contrast to Galactic black hole 
binaries, many ULXs have a cool (~ 0.1 — 0.2 keV) ac- 
cretion disk suggesting a black hole mass of > 100A/ Q 
(e.g. Kaaret et al. 2003; Miller, Fabian & Miller 2004a; 
Winter, Mushotzky & Reynolds 2006). While Galactic 
black hole binaries have a good correlation between the 
spectral shape and luminosity, the behaviors of ULXs 
are more complex (e.g. Winter et al. 2006; Kong et 
al. 2007; Feng & Kaaret 2009; Kajava & Poutanen 
2009; Vierdayanti et al. 2010). ULXs seem to require 
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more complicated accretion geometry invoking outflows, 
corona, massive donors, and super-Eddington accretion 
flows (e.g. Stobbart et al. 2006; Poutanen et al. 2007; 
Patruno & Zampieri 2008; Gladstone et al. 2009). 

It is possible that ULX is a distinct class of systems 
comparing to Galactic black hole binaries. Apart from 
a few very bright ULXs (also known as hyperluminous 
X-ray sources, e.g. Farrell et al. 2009), a > 100M© 
black hole may not require based on the observed lu- 
minosity. Instead, black holes with a few tens of solar 
masses may be more likely. Theoretical models involv- 
ing binary mergers (Belczynski et al. 2004) and low- 
metallicity massive (~ 40 — 5OM ) progenitors (Mapelli 
et al. 2009; Zampieri & Roberts 2009) are proposed to 
explain a population of ULXs with ~ 30 — 9OM black 
holes. 

Alternatively, ULX could be a typical stellar-mass 
black hole with geometrically or relativistically beamed 
emission (King et al. 2001; Kording et al. 2002) so 
that the X-ray luminosity does not exceed the Edding- 
ton limit. Furthermore, the stellar-mass black hole may 
in fact accrete materials at or above Eddington limit via 
a slim disk (Ebisawa et al. 2003), or a radiation pressure- 
dominated accretion disk model (Begelman 2002). It is 
also possible due to combination of both scenarios (King 
2008). 

Holmberg IX X-l is a famous ULX located near the 
galaxy M81 and it is about 2 arcmin from M81's dwarf 
companion, Holmberg IX. The source was first discovered 
by the Einstein Observatory (Fabbiano 1988) and has 
been observed by all major X-ray observatories through- 
out the last 20 years (La Parola et al. 2001). Apart from 
the X-ray flux variability, Holmberg IX X-l is also one 
of the first ULXs shown to have a cool (~ 0.1 — 0.2 keV) 
accretion disk, leading to a suggestion of an intermediate- 
mass black hole accretor (Miller et al. 2004a). It is 
proven that monitoring observations can reveal the phys- 
ical nature of Galactic X-ray binaries by tracking their 
flux and spectral evolution as well as their correlation 
(see Remillard & McClintock 2006). Until now, it has 
been quite difficult to monitor ULXs due to their dis- 
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tances and crowding location. Swift is the first X-ray 
telescope with reasonable spatial resolution and sensitiv- 
ity to perform such observations. Recently, Kaaret & 
Feng (2009) reported Swift monitoring observations of 
several ULXs including Holmberg IX X-l, NGC 5408 X- 
1, and NGC 4395 X-2. 

In this paper, we report a Swift monitoring observation 
of the ULX, Holmberg IX X-l with a focus on the spec- 
tral behaviors. We describe the observations and data 
reduction method in §2. The results are present in §3, 
and a discussion is in §4. 

2. OBSERVATIONS AND DATA REDUCTION 

Holmberg IX X-l is one of the ULXs monitored with 
the X-ray Telescope (XRT) of Swift (Kaaret & Feng 
2009). The source has been observed with Swift regu- 
larly since 2006. In particular, there is a guest observ- 
ing program (PI: Kaaret) for an intensive monitoring of 
Holmberg IX X-l since 2008 December. Furthermore, we 
proposed a follow-up monitoring program in mid- April 
2009 and all the data obtained between April 24, 2009 
and 2009 July 23 are from this new program but with 
the same ObsID (90008). In this work, we focus on the 
data obtained after 2008 December. In addition, we also 
include data from ObsID 90079 obtained between 2009 
April and 2010 March. We only used data taken in pho- 
ton counting mode. Timing analysis and hardness vari- 
ation for data taken between 2006 July and 2009 June 
have been reported in Kaaret & Feng (2009). During the 
period we are interested, we obtained 177 XRT observa- 
tions of Holmberg IX X-l with a total exposure time of 
242.1 ks. 

We extracted XRT light curves and spectra of Holm- 
berg IX X-l by using the XRT products generator., 
(Evans et al. 2007, 2009). In brief, this software first 
creates an image from the event list and identifies our 
target for each observation. Only events with energy in 
the range 0.3-10 keV with grades 0-12 are included. A 
circular source extraction region is chosen to maximize 
the signal-to-noise ratio of the light curve and spectrum. 
For the background, an annulus centered on the source is 
used. Then for each observation, source and background 
counts are extracted. Source count rates are corrected 
for the good time interval, and losses due to bad pixels 
and bad columns. The source spectrum is extracted from 
a combined event list of all the observations considered. 
Similarly, the software automatically extracts and com- 
bines all the ancillary response files (ARFs) ; each ARF is 
weighted according to the proportion of counts in the to- 
tal source spectrum. Background spectrum is extracted 
from an annulus region centered on the target exclud- 
ing any sources in the extraction region. An appropriate 
redistribution matrix file is selected from the CALDB. 

We grouped the source spectrum with at least 20 
counts per spectral bin before fitting in order to allow % 2 
statistics for finding the best-fitting parameters. We per- 
formed spectral analysis using the HEAsoft X-ray spec- 
tral fitting package XSPEC version 12.5. All spectral fits 
were made in the 0.3-10.0 keV band. All errors reported 
in this work are 90% confidence errors. 

3. ANALYSIS AND RESULTS 
6 http:/ /www. swift. ac.uk/user_objccts 



The X-ray long-term light curve of Holmberg IX X-l in 
the 0.3-10 keV is shown in Fig. 1. Our light curve created 
by the XRT products generator is almost identical to the 
one generated by Kaaret & Feng (2009), except that our 
count rate is systematically slightly higher than theirs. It 
is likely due to different energy range as well as different 
correction methods for bad pixels and the point-spread- 
function. It is very clear from Fig. 1 that the source 
exhibits substantial variability. In order to search for 
any modulation in the light curves, we used the Lomb- 
Scargle periodogram (LSP; Lomb 1976; Scargle 1982), 
a modification of the discrete Fourier transform which is 
generalized to the case of uneven spacing. Like Kaaret & 
Feng (2009), we found two peaks at ~ 20 days and ~ 60 
days in the LSP. However, both peaks are not statistically 
significant with a maximum power of 6.6 near a period 
of 19 days (the 99.9% significance level has a power of 
11.9). The light curve also shows three distinct intensity 
states. From 2008 December to 2009 early- April, the 
source count rate varies between 0.15 counts per sec and 
0.3 counts per sec, and shows obvious modulation. The 
source becomes fainter from 2009 mid-April (~ Day 120 
in Fig. 1) to 2010 early- January with an average count 
rate of 0.1 counts per sec; we defined this period as the 
"low" state. Since 2010 January 10, the source intensity 
has increased to 0.3-0.4 counts per sec and we name it 
as the "high" state. 

Because of the low count rate for each Swift observa- 
tion, we added all the spectra in similar state together to 
study the spectral behavior. We divided the light curve 
into three parts: 1) the "variable" state with data taken 
between 2008 December 1 and 2009 April 1; 2) the "low" 
state with data taken between 2009 April 23 and 2010 
January 3, and 3) the "high" state with data taken be- 
tween 2010 January 10 and 2010 March 31. We also con- 
sidered a co-added spectrum from all the data. We first 
fitted all the spectra with an absorbed power-law model 
and the spectral parameters are listed in Table 1. All 
spectra cannot be fitted satisfactorily with a power-law 
model. Similarly, an absorbed MCD blackbody model is 
not an acceptable fit to any of the spectra when it is the 
only continuum component. 

We next considered to apply a MCD plus power-law 
model that provides good fits to a sample of ULXs (e.g. 
Kaaret et al. 2003; Miller, Fabian k Miller 2004a; Win- 
ter, Mushotzky & Reynolds 2006). The additional MCD 
component is statistically significant for all states. For 
the "low" state spectrum, the addition of a disk com- 
ponent is significant at the > 4c level of confidence. In 
particular, the best-fitting disk temperature is very low 
with kT — 0.19 keV (see Table 1), which is consistent 
with previous observations (Miller et al. 2004b) . For the 
"variable" and "high" state spectra, although the addi- 
tional MCD component is statistically required, the spec- 
tral parameters are completely different comparing to the 
"low" state spectrum. The disk temperature is much 
higher (1.69-2.25 keV against 0.19 keV) and the photon 
index is also very different (2.65-3.47 against 1.68). Fur- 
thermore, this model is only marginally acceptable for 
the variable states. 

We also fit a dual thermal model consisting of a cool 
blackbody continuum at low energies and a hot disk 
blackbody component at high energies (Stobbart et al. 
2006). This model provides the best fitting for all states 
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Fig. 1. 



Swift XRT 0.3-10 kcV lightcurvc of Holmbcrg IX X-l between 2008 December 1 and 2010 March 31. 



as well as all the data combined with very similar spec- 
tral parameters (kT = 0.2 keV and kT in — 2 keV; Table 
1 and Figure 2). 

Finally we consider a more physical self-consistent 
Comptonization spectrum using the DISKPN+EQPAIR 
model (Stobbart et al. 2006; Gladstone ct al. 2009). 
The EQPAIR model (Coppi 1999) allows thermal and 
non-thermal electron distributions. We tie the temper- 
ature of the seed photons to that of the inner accretion 
disk described by the DISKPN model. This model can 
describe the "high" state data equally well as the dual 
thermal model with a disk temperature similar to that 
of the cool blackbody component of the dual thermal 
model. For the "low" and "variable" states, the fits are 
slightly worse than that of the dual thermal model. 

4. DISCUSSION 

We obtained a long-term X-ray light curve of Holmberg 
IX X-l by using Swift XRT and found that the source 
transited from a "variable" state to a "low" state, and 
then back to a "high" state. The co-added spectra of 
all states can be marginally described with a MCD plus 
power-law model, and a dual thermal model (blackbody 
plus MCD) provides the best fits. For the MCD plus 
power-law model, the "low" state spectrum can be best 
fitted with a cool accretion disk (kT in = 0.19 keV) plus a 
power-law (r = 1.7). However, during the "variable" and 
"high" states, the spectral parameters are completely dif- 
ferent (see Table 1). This is very unusual given that 
the flux difference is only about a factor of 2-3. When 
we examined the hardness ratios, we did not find any 



significant change throughout the whole monitoring ob- 
servation. Kaaret and Feng (2009) also showed by using 
hardness ratios that the source did not exhibit significant 
spectral change. We therefore conclude that the appar- 
ent spectral change based on spectral fitting may not 
be real, and it is simply a swap between the MCD and 
power-law components. More specifically, one common 
ULX model consists of a cool MCD and a hard power- 
law just like the spectrum in the "low" state; during 
the "variable" and "high" states, we have a soft (steep) 
power-law and a hard MCD. In fact, when we fixed the 
MCD temperature of the "variable" and "high" state 
spectra at the "low" state value (0.19 keV), the other 
spectral parameters are consistent with the "low" state. 
Moreover, Stobbart et al. (2006) rejected the soft power- 
law plus hard MCD model using XMM-Newton data. 
Thus we do not consider there is a spectral change during 
the Swift monitoring campaign. 

While a MCD plus power-law model does not give a 
consistent result in all the states, the dual thermal model 
fits the data quite well and it is also the best model 
among all the spectral models we considered (see Table 
1 and Figure 2). This dual thermal model is motivated 
by the presence of optically thick outflowing winds from 
a black hole accreting at or above the Eddington limit 
(King & Pounds 2003). When this happens, the thick 
wind develops a photosphere emitting at a certain tem- 
perature. This blackbody temperature may explain the 
ultrasoft (~ 0.1 keV) X-ray component of some ULXs. 
More recently, Stobbart et al. (2006) applied this dual 
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TABLE 1 
Spectral fit parameters 



Model parameter 


Low state 


Variable state 


High state 


All data 


Power-law model 


N H (10 2i cm" 2 ) 


1.74 ±0.13 


2.16 ±0.10 


2.96 ±0.23 


2.16 ±0.07 


r 


1.74 ±0.04 


1.86 ±0.03 


2.00 ± 0.06 


1.86 ±0.02 


X 2 /dof 


430.6/315 


537.5/445 


238.2/220 


786.7/535 


io.5-10 (lO^ergs" 1 ) 


1.00 


1.66 


2.48 


1.58 


MCD + Power-law model 


N H (10 21 cm" 2 ) 


2 60 +u ' ti4 
Z - DU -0.51 


2.80t£E 


4.12±i;«{ 


°' lu -0.44 


r 

kT m (keV) 
Norm 


1.68 ±0.06 

20+ 06 
u - zu -0.05 
71 +206 
— 55 


2 65+ - 53 

9 9t-+0.23 

Z - ZJ -0.18 

012+ 0007 
U ' UAZ — 0.006 


3.47 ± 1.03 

1 69+ 018 

1 - oy -0.15 

061 + 038 
U.UD1_ 027 


3 02+ ' 41 
9 1 1 +0.12 

z - 1J --o.ii 

01 6+ 005 
U.U1D_ 004 


X 2 /dof 

FMCD / FTotai a 


410.8/313 
0.13 


515.7/443 
0.46 


205.7/218 
0.50 


695.4/533 
0.52 


Lo.5-10 (lO^'crgs- 1 ) 


1.15 


1.75 


3.02 


1.78 



MCD + Blackbody model 



N H (10 21 cm- 2 ) 


1 39+ U ' iB 
1 - oa -0.18 


1.20 ±0.10 


l.OD_ 2g 


1 25+ Uil 


kT (keV) 


n 99+O.O2 
u - zz -0.12 


0.26 ±0.01 


0.26 ± 0.04 


0.25 ±0.01 


kT m (keV) 


2 05+ ' 12 

z - uo -o.io 


1.86 ±0.05 


1 60+ 010 
1 - DU -0.03 


1.84 ±0.06 


Norm 


014+ - 001 
U-U — 0.003 


0.037 ± 0.004 


086+ ' 028 
U.UOD_ 023 


0.033 ± 0.004 


X 2 /dof 


348.2/313 


489.7/443 


203.7/218 


625.3/533 


F M CD/F Tot ai a 


0.85 


0.85 


0.86 


0.85 


£0.5-10 (lO^'crgs- 1 ) 


0.88 


1.38 


2.0 


1.35 



DISKPN + EQPAIR model 

0+0.21 

-0.09 



9 9f? + W-3» 
Z - ZO -0.36 
,+0.11 
-0.06 
-+1.06 
-0.88 



, +U.1U 

-0.23 



N H (10 21 cm- 2 ) 
T max (keV) 
lh/h b 

T 

X 2 /dof 
£0.5-10 (10 40 ergs' 1 ) 



2.22 ±0.25 
1-0.06 
-0.01 



0.26Z 



cr 41 0.51 
°- 41 -0.25 

26.66li;°| 
354.8/312 
1.03 



1.92] 
0.28 ±0.02 
3.27 ±0.13 
19.45 ± 1.18 

496.7/442 
1.72 



2.01: 

0.31+0.06 0.29 ±0.02 

3-05-^0.88 

23.74 ±6.24 23.77! 
203.9/217 643.9/532 
2.17 1.47 



a +o.io 

="-0.14 
7+0.91 
-0.99 



NOTES. — All quoted errors are 90% confidence. Luminosities arc calculated assuming a distance of 3.6 Mpc. 
a Flux ratio between MCD component and total flux. 

b Ratio between the compactness of electron and the compactness of the seed photon distribution. 



thermal model to describe the XMM-Newton spectra of a 
sample of ULXs (Holmberg IX X-l is one of the sources) 
and found out that 10 out of 13 ULXs can be fitted satis- 
factorily, with kT = 0.15 - 0.3 keV and kT in = 0.8 - 2.2 
keV. This indicates that we may be observing both the 
accretion disk as well as the wind from the central black 
hole. In this scenario, these ULXs are simply the exten- 
sion of stellar-mass black hole accreting at or above the 
Eddington limit. 

Lastly, the "high" state spectrum can also be fit with 
a cool accretion disk plus Comptonized corona model. 
The best-fit parameters are consistent with other ULXs 
observed with XMM-Newton (Gladstone et al. 2009). 
Based on the fitting statistics (see Table 1), we cannot 
distinguish the dual thermal model and accretion disk 
plus Comptonized corona model. For the "low" and 
"variable" states, the accretion disk plus Comptonized 
corona model is slightly worse than the dual thermal 
model and we cannot totally rule out this model. It 
is therefore not clear if the source underwent a spectral 
change. Using similar approach and employing a disk 
plus thermal Comptonization model, Vierdayanti ct al. 
(2010) assert that spectrum evolves with a definite trend 
in which the corona temperature decreases and its opti- 
cal depth increases as the source becomes brighter. Our 
data do not support this conclusion. Indeed, Vierdayanti 
ct al. (2010) present no statistical tests of the correla- 
tion of their fitted spectral parameters with luminosity 
to support this conclusion. However, there may be some 



subtle spectral differences in the dual thermal model. For 
instance, while the ultrasoft soft X-ray component is sim- 
ilar, the best-fit kTi n of the "high" state is slightly lower 
than that of the other two states. 

Holmberg IX X-l has been observed with XMM- 
Newton several times and indeed it is one of the first 
ULXs to test the MCD plus power-law model (Miller 
ct al. 2004b). In Stobbart et al. (2006), the MCD 
plus power-law model provides a better fit than the dual 
thermal model although both are statistically acceptable. 
In our Swift monitoring observation, the dual thermal 
model (and perhaps also the accretion disk plus Comp- 
tonized corona model) can always provide a better fit 
while the MCD plus power-law model is only acceptable 
in the "high" state data. If this is true, Holmberg IX X-l 
may be a black hole accreting at or above the Edding- 
ton limit, instead of a massive black hole with > 100M Q . 
This model is also supported by a study of Holmberg IX 
X-l using XMM-Newton and ASCA data in which the 
spectra can be described by a slim disk model indicat- 
ing that the source is accreting near the Eddington limit 
(Tsunoda et al. 2006). 

In our dual thermal model, the disk component con- 
tributes about 80% of the total emission which is totally 
different comparing to the cool disk spectra (in which 
the MCD component contributes only ~ 20%) discussed 
by Miller et al. (2004b). Therefore it is likely that 
the source during our Swift monitoring as well as the 
XMM-Newton observation taken in 2001 and ASCA ob- 
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Fig. 2. — Swift XRT 0.3—10 keV spectra of Holmbcrg IX X-l taken during the "low" state (upper left), "variable" state (upper right), 
"high" state (lower left), and all data (lower right). All spectra can be well described by a dual thermal model with a blackbody temperature 
of ~ 0.2 keV and a MCD temperature of ~ 2 keV (See Table 1 for spectral parameters). 



servation taken in 1999 (Tsunoda et al. 2006) is in a 
disk-dominant state. We re-analyzed the XMM-Newton 
data used in Tsunoda et al. (2006) and found that the 
spectrum can also be well described by the dual thermal 
model with kT = 0.32 keV and kT m = 1.56 keV; this 
is consistent with the "high" state of our Swift obser- 
vations. We also modeled the XMM-Newton data with 
the DISKPN+EQPAIR model and the best-fit parame- 
ters {T max ~ 0.22 keV; r ~ 20) are similar to those of 
the Swift "high" state. 

The observed spectral parameters of this dual ther- 
mal model are roughly consistent with the calculation 
of Poutancn et al. (2007) in which the outflow (low 
temperature component) is from the spherization ra- 
dius. This also applies to the cool disk component of the 
DISKPN+EQPAIR model. Following Poutanen et al. 
(2007), a blackbody temperature of 0.2 keV corresponds 
to an accretion rate of about 16 times of the Eddington 
rate for a 100M Q black hole. For the high temperature 
component from the inner disk, the predicted tempera- 
ture is about 1 keV but Poutanen et al. (2007) noted 
that it can be up to 4 keV for a rotating black hole. For 
Holmberg IX X-l, the disk temperature is about 2 keV. 
Following Feng & Kaaret (2005), the fractional Edding- 
ton luminosity can be estimated as 



kT 



L 



x 



1.2keV/ \1.3 x 10 39 ergs s" 



0.5 



(1) 



where f3 is determined by the black hole spin with 
8 = 1 for a Schwarzschild black hole and 1 /6 for a max- 
imally rotating Kerr black hole. Given an X-ray lumi- 



nosity of <~ 10 40 ergs s i , Holmberg IX X-l appears to 



be above the Eddington limit by a factor of 1.3 even for 
a maximally rotating black hole. If the source radiates 
at the Eddington limit, we estimate that the mass of the 
black hole is about 100M Q based on the observed lumi- 
nosity. If Holmberg IX X-l is a non-rotating black hole, 
the fractional Eddington luminosity can be as large as 
7 and the disk temperature corresponds to a ~ 10M Q 
black hole. It is therefore reasonable to expect that the 
black hole mass of Holmberg IX X-l is between 10 and 
100 M Q . It is worth noting that a - 30 - 90M Q black 
hole can be formed via binary mergers (Belczynski et 
al. 2004) or a massive progenitor in a low metallicity 
environment (Mapelli et al. 2009; Zampieri & Roberts 
2009). 

However, as noted by Stobbart et al. (2006), this dual 
thermal model requires a specific geometry and viewing 
angle so that both an optically thick outflow and the 
inner regions of the accretion disk can be seen at the 
same time. From recent optical observations, Holmberg 
IX X-l is surrounded by a shell- like shock ionized nebula 
(Miller 1995; Wang 2002; Pakull & Grise 2008) and this 
is consistent with our interpretation of the X-ray spectra 
that strong wind puffs from the disk. 

In summary, the nature of the compact object of 
Holmbcrg IX X-l is still a mystery. Previous XMM- 
Newton observations suggest that it is an intermediate- 
mass black hole based on a cool accretion disk model. 
Timing study shows that it may be a 50 — 200M Q black 
hole (Dewangan et al. 2006). Our study using the long- 
term Swift monitoring observations indicates that a dual 
thermal model can always provide the best fit suggest- 
ing a 10 — 100M Q black hole accreting at or above the 
Eddington limit. In order to investigate if there is any 
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spectral change and also the nature of the black hole, 
Holmberg IX X-l is therefore deserved for further study 
using monitoring campaign. 
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